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I. INTRODUCTION 


For quite some time, physicists have considered Einstein’s General Relativity (GR) to 
be an effective theory of gravity. Therefore, in order to find the happy marriage between 
quantum theory and gravity, we need to know the underlying fundamental theory of gravity. 
One recent proposal in this address has been the Horava-Lifshitz (HL) theory [l|, that is, a 
power-counting renormalizable theory with consistent ultraviolet (UV) behaviour. Further¬ 
more, the theory has one fixed point in the infrared (IR) limit namely GR [1-4 ]. In terms of 
the above, black holes (BH) are important solutions for field equations in any gravitational 
theory, including those such as Einstein-Hilbert, Brans-Dicke, HL, f(R), string theories and 
any generalisation or modification of Einstein’s gravity. At the quantum level, BH play the 
same role as hydrogen atom and we hope they give us some clues about the observables 
of any quantum theory of gravity. As such, it is important to study the physical prop¬ 
erties of BH solutions, such as decay rate, greybody factors, or their quasinormal modes. 
Quasinormal modes (QNMs), known as “ringing” in BH, are very important in order to un¬ 
derstand the classical and quantum aspects of BH physics. The QNMs give us hints about 
the stability of BH under consideration, as in this paper, and can be used to compute the 
spectrum of the area operator using the semiclassical approach developed by Hod jfj . The 
determination of QNMs is based on the dynamics of matter fields and on the metric pertur¬ 
bations in the BH background. In this work indeed, we are interested in the stability of the 
1 + 1-dilatonic HL BH using a QNMs’ approach; QNMs associated with the perturbations of 


different fie’ 
AdS space 


7- 


ds have been considered in different works jf]], including those involving dS and 


14] and higher dimensional models, where the QNMs can be computed for a 


brane situated in the vicinity of a D -dimensional BH [15j]. A similar situation occurs in 2 +1 


dimensions 


16H18|. and for acoustic BH 19r21[. QNMs of dilatonic BH in 3 + 1 dimensions 


can be found in Refs 
much attention 


22 


241 ]. Two-dimensional theories of gravity have recently attracted 


25rl27l] as simple toy models that possess many features of gravity in higher 


dimensions. They also have BH solutions which play important roles in revealing various 
aspects of spacetime geometry and quantization of gravity, and are also related to string 


studied in 


30] 


theory 28, 2J- The QNMs of 1 +1 dilatonic BH for scalar and fermionic perturbations were 


33]. 


The determination of QNMs for a specific geometry implies solving the field equations for 
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different types of perturbations (scalar, fermionic, vectorial, etc.), with suitable boundary 
conditions that reflect the fact that this geometry describes a BH. The QNMs of a classical 
scalar perturbation of BH are defined as the solutions of the Klcin-Gordon equation char¬ 
acterised by purely ingoing waves at the horizon, <f> ~ e -uu(t+r)^ s j ncej least classically, an 
outgoing flux is not allowed at the horizon. In addition, one has to impose boundary con¬ 
ditions on the solutions in the asymptotic region (infinity), and for that reason it is crucial 
to use asymptotic geometry for the spacetime under study. In the case of an asymptotically 
flat spacetime, the condition we need to impose over the wave function is to have a purely 
outgoing wave function $ ~ e -iw(t-r) a £ ^he i n fi n ity [7], In general, the QNMs are given 
by ujqnm = Ur + iojj, where u>r and up- are the real and imaginary parts of the frequency 
uqnm, respectively. Therefore, the study of QNMs can be implemented as one possible 
simple alternative test for studying the stability of the system. In this sense, any imaginary 
frequency with the wrong sign would mean an exponentially growing mode, rather than a 
damping one. 

The organisation of this article is as follows: In Section II, we describe briefly the HL 
theory and specify the 1 + 1-dilatonic BH solutions. In Section III, we compute the QNMs 
and explore the criteria for the stability of the two BH metrics under consideration. We 
finish with conclusions in Section IV. 


II. GENERALITIES OF THE HORAVA-LIFSHITZ GRAVITY 

In the following, we will describe the HL theory as developed in Ref. [41]. The HL theory 
provided a new approach to quantum gravity and its principal idea is based on the breaking 
of the Lorentz invariance by equipping the spacetime with additional geometric structure, a 
prefered foliation which defines the splitting of the coordinates into space and time; in this 
theory the Lorentz invariance is assumed to appear only at the low energies limit. One can 
decompose the spacetime as follows 


ds 2 = ( N 2 - NiN^dt 2 - 2 Nidx'dt - h ij dx i dx j , 


( 1 ) 


where N , N l are the lapse and shift functions respectively and is the three-dimensional 
metric. The action is given by 


s=Mk 


d 3 xdtVhN ( KijK ij - 7 K 2 - V) , 


( 2 ) 
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where Mpi is the Planck mass, 7 is a dimensionless constant and K i3 is the well-known 
extrinsic curvature tensor, which is stated in the ADM formulation as 

Kij = (hij — VjA? — VjAj^ , (3) 

being K its trace. The last term in (l2j) . V, is invariant under three-dimensional diffeo- 
morphisms and is known as the “potential” term. This term is a function of the three- 
dimensional metric and its derivatives. In explicit form we have 

V = -^ R+ Jp ^ lAR + ^ Rij + -) + W- {°i a2r + ^R,,R 3k R\ + ...)> ( 4 ) 

where £, 7 r n , a n are coupling constants, Rij and R are the Ricci tensor and the scalar curvature 
constructed with the spatial metric. A := hAViVy The introduction of the “potential” 
term in (J2]) improve the UV behavior of the graviton propagator and additionally leads to 
different scaling of space and time 

x —^ p 1 x, t —^ p A —^ A, A* —^ p^Aj, hij —^ h^. (5) 


When the lapse function depend only on time, A = N(t), we say that we are dealing with 
the “projectable” version of the HL theory and the “non-projectable” version is given when 


the lapse function can depend on space and time. In Ref. 


34J an extension of the non- 


projectable version of HL gravity was made by the introduction of a extra mode in the 
“potential” term, i.e., V(hij ) —> It was shown that this extra mode can acquire a 

regular quadratic Lagrangian. The extra mode is given by 

d,N 




N 


( 6 ) 


Geometrically this vector represents the proper acceleration of the unit normals to the 
spatial slices. For the two-dimensional case there are only two terms that contribute to the 
quadratic Lagrangian: R and aia\ 


A. Lowest dimensional Horava-Lifshitz Black Hole 

The HL gravity has two-dimensional solutions that characterise dilatonic BH and can be 
used to study the physical properties of BH in general; furthermore, some features of this 
theory, due to the fact that it utilizes two dimensions, open the possibility of understanding 
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physical consequences in higher dimensional theories. As a summary, let us start with the 


HL-dilaton gravity in two dimensions presented in Ref. 


35), 


S = Shl + Ss, 


(7) 


where, as mentioned before, the quadratic Lagrangian for the HL theory in two dimensions 
comes from the contribution of the terms R and a*a* 

Shl = ^ 7 p J dtdxy/g ((1 - X)K 2 + rjg u a icq) , (8) 

and 


Sep — I dtdxNy/g 


— (dt<t> - iWi) 2 - «(Vx0) 2 - V{<j>) - /30VV - 10 


, (9) 


where a, /3, rj and ? are constants. Using the fact that K = 0 and admitting Ni = 0 together 
with the relativistic limit f3 = s = 0, we are left with the action 




dtdx ( --rjN 2 a\ + aN 2 (\j 2 — V(<f>) 


( 10 ) 


From now on, the prime denotes derivative with respect the coordinate x. In two dimensions 
the extra mode a, is simply a\ = d\ In N = (In N) . In Ref. 35], a new set of BH solutions 
in two-dimensional HL gravity was found for action (ITU]) . The solutions are described by 


N(x) — 4 /— x 2 — 2C\X H-h ——— + 2 C 2 , 

r] Tjx 6r]x z 


( 11 ) 


and 


4>{x) = bi \ / —x 2 - 2Cix + — + 7 -^— + 2 C 2 , 

7] rjx 61JX Z 


( 12 ) 


these solutions were obtained by using the quantity 

v *( x ) = A + 73 + 73- 


(13) 


the derivative of the scalar potential given as a function of an implicit scalar field which in 
turns depends on the spatial coordinate. This was done because for generalized potentials 
is not always possible obtain analytical solutions. We would like to focus our attention on 
the following two cases, 
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• First case: Described by fixing the constants in the following way; A — B — C — 0, 
C\ = —M, C ‘2 = —1/2 and 77 = 1. In this case = 0. Therefore, the metric for this solution 
can be written as follows 


ds 2 = —(2 Mx — 1 )dt 2 + 


- dx 2 , 


(14) 


2 Mx - 1 

where the parameter M is related to the lapse function N. This solution was found for the 
first time in 


• Second Case: Here we fix the constants as A = A, B = C = 0, C\ = — M and C 2 = — §; 
therefore, we have — A, and the solution is given by 


ds~ = 


— 1 x 1 + 2 Mx 


l\V 


dt 2 + 


■0 x 2 + 2Mx 


-dx 2 . 


(15) 


The horizon of the black hole is located at 

r/M 


x± = 


A 


± 


/ 77 / rjM 2 ^ 


If we define the variables u = \J A /-qx + \Jr\/KM and = V(r//A)M 2 + e one gets 

l 2 


ds 2 = — (u 2 — u \) dt 2 + 


-du , 


(16) 


(17) 


(u 2 - u 2 ) 

which is a suitable expression to study the quasinormal modes of this black hole, and we 
have defined l = (A/ 77 ) 1 / 4 . I n this new coordinate system, the horizon of the black ho’ 
located at u — u + . The spacetimes, described by (TT4|) and (TlST) . are conformally flat 37]. 


e is 


III. QUASINORMAL MODES 

In order to study the QNMs, we consider a scalar field minimally coupled to gravity 
propagating in the background of the two-dimensional HL BH. We consider the following 
action for the scalar field 

S[<p] = I d 2 x^g ^-^(Vy ?) 2 - -m 2 (p 2 ^j , (18) 

where m is the mass of the scalar field. From the variation of 5(J) the field equation is given 
by 

□79 — m 2 (p = 0, (19) 

where □ is the D’Alambertian operator, in the following sections, we will solve the Klein- 
Gordon equation dT9]) for the spacetimes described in the previous section. 
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A. Spacetime metric ds 2 = —(2Mx — 1 )dt 2 + 


1 

2Mx—l 


dx 2 


1. Massive Scalar Field 


The case of a massive scalar field perturbing the background described by the metric (THj) 


was discussed for first time in 
(USD, is represented by 


36|,l37]. If we use (p(t,x) = e lut ip(x) the equation of motion 


(2 Mx -1 )FM. + + 


u 


dx 2 


dx 


2 Mx - 1 


ip(x) — m 2 (p(x) = 0. 


( 20 ) 


For this metric and massive scalar field, the QNMs were computed in 38|, where the authors 
claimed that the solution for the QNMs are completely different from the standard case where 
the QNMs have a discrete spectrum. They proposed a real and continuous spectrum for the 
QNMs of a scalar perturbation. In the next section, we look at this situation, but instead 
use the confluent hypergeometric function 0 Fi(a, 6; x) place of the modified Bessel functions. 
Now, using the tortoise coordinate defined by x* = ^ ln(2Mx — 1) we can write the Eq. 
(|25D as a Schroedinger type equation with effective potential 1 

V ef f = m 2 e 2Mx \ 


( 21 ) 


this potential diverges when x* —> oo. If we consider the variable z = jje Mx * 36|, after some 


straightforward algebra, the equation of motion 


MMlL + zMd + + ,2) v{z) = o, 


can be written as the Bessel equation, 

( 22 ) 


dz 2 ' dz 

where v — this equation can be transformed into the confluent hypergeometric equation 
using the change (p(z) = (2 iz) v e~ 2z F(z) 39], and we find 


+< 2v + 1 - + 2i ("+ \ I F ^ - °’ 


(23) 


whose solution is given in terms of the confluent hypergeometric, or Kummer, functions 


1 1 
F{z) = A$(v + -, 2 u + 1; 2 iz) + B(2iz)~ 2v ^(-v + -, 1 - 2z/; 2 iz), 


(24) 


where A, B are constants and the number v in general need not be an integer 39]. In the 


following we will consider two cases for the v parameter since we are interested in exploring 
all its possible values. 
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• non-integer v 

In order to compute the QNMs, we need to impose adequate boundary conditions that 
represent a purely outgoing wave at infinity and purely ingoing wave near the horizon of the 
BH (a condition often used in flat spacetime). Another situation occurs when the asymptotic 
behaviour of the spacetime is not flat, e.g. asymptotically AdS space; in these kind of spaces, 
the potential diverges at infinity, and we can therefore impose p = 0 (Dirichlet boundary 
condition) or ^ = 0 (Neumann boundary condition) at infinity. As we can see from Eq. fj2Tjl . 
in this case we have an asymptotically AdS space. Therefore, we need to apply boundary 
conditions to the QNMs over the general solution of Eq. (J23]), which is given by 

p(z) = Ae~ 2z z v <h(z/ + i, 2z/ + 1; 2 iz) + Be^ 2z z~ u ^(—v + 1 — 2z/; 2 iz), (25) 


to satisfy the boundary conditions properly, we set A = 0 in order to have only ingoing 
waves at the horizon (2 = 0 ). The asymptotic behavior of <f> at infinity is given by 40] 


r (p) 


(26) 


and therefore our solution at the infinity (z —>■ 00 ) reads as follow, 

p(z) ~ Bz~^e A 


- i e 2(i-i) z r (! 2z/ ) 

T(—v + 1/2)' 


(27) 


We can see that the scalar field vanishes as z —> 00 , this confirms the absence of QNMs 
for this HL BH under scalar perturbations, a similar situation was found in Ref. 18] for 


the QNMs of the extremal BTZ BH. The conclusion of this case was discussed in 


and represents a continuous spectrum. Additionaly, if we impose the Neumann boundary 
condition, we obtain a similar asymptotic vanishing behavior for the flux 


t/n */ sd<p(z) ( ^dip*(z) 
Jr{z) = P {Z )—7- P{Z)- 


(28) 


dz dz 

In light of the meaning of QNMs, for any BH perturbation its geometry produces damped 
oscillations, this is the so-called ringing in BH. It is well known that the frequencies of 
these oscillations and their damping periods are completely fixed by the BH properties, 


and as such, are independent of the nature of the initial perturbation. In Ref. 


it was 


shown that, for scalar perturbations the oscillations have a continuous spectrum and are 
not discrete, as would be expected for a BH. This result is very strange and in our opinion, 
devoid of physical meaning; this is because it is well known that oscillations of QNMs are 
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similar to normal modes of a closed system. In the next section, we consider the second 
solution of the confluent hypergeometric equation and show that QNM oscillations have a 
discrete spectrum. 

• v integer 

Now we are considering the case where v is an integer number. For v integer, the solution 
of the Eq. (1231) changes. In the expression (l25]h the function <F(p, q;z) of the second term 


must be replaced by 


41] 


71=1 


where 


r (n + a)T( r y)B n z r ‘ 

' r(a)r(n + 7)n! 

r( 7 )r?r + a — 7 + ir (7 — n — 1)(—l) n 
T (a)n\z ' r ~ n ~ 1 

Via) 


W(a, i;z) — M(a,xz) (in z + ip( 1 - a) - ^( 7 ) + C) + ^ 

OO 

+ (-D’E 


+ 


i>(a) = 


F(a)’ 


represents the digamma function, C = 0.577216... is the Euler’s constant, and 

1 \ (\ 1 1 


B n = 

Then, we have 


1 1 

—I-— + ... + 

a a + 1 


a + n — 1 


1 1 

—I-— + ... + 

7 7 + 1 


7 + n — 1 


(p(z) = Ae lz z v ^{v + -, 2v + 1; 2iz) + Be lz z v W(p,q;z). 


(29) 


(30) 


(31) 


In order to have only ingoing waves at the horizon (z = 0), we set A = 0. The asymptotic 
behaviour of the IF-function at infinity is given by 


W(p,q;y ) -+ 7 rcot(vrp)Q^?/ p q e y . 

r (p) 

Therefore, the general solution at inhnity reads as follows, 

1 


<p{z) ~ Bn cot 


—7T\U + 


(41/+I) T(2z/+1) 2(i-l)z 

r(-i/-1/2) 


(32) 


(33) 


If we consider that u + 1/2 = 2n + 1 , where n is an integer number, we are able to fulfill the 
Dirichlet boundary condition at inhnity. From this result we can obtain the frequency of 
the QNMs as 

oj = —inM, (34) 

and using the Neumann condition for a vanishing flux at inhnity, we obtain the same result 
for the QNMs as expressed in Eq. 
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2. Massless Case 


For the metric CUD, when m = 0, we use the standard definition for QNMs, the Klein- 
Gordon equation (flUj) which reads 


(2 Mx 


dV(x) 

1 dx 2 


+ 2 M 


dip(x) 

dx 


+ 


oj 


2 Mx - 1 


<f(x) = 0, 


(35) 


where we have assumed <p(t,x) = <p(x)e lwt . If we define the quantity x + = 1/2 M and the 
change of variable z — 1 — x + /x we can write equation (1351) as follows 




d 2 ip 
dz 2 




(1 — z) dip 

- - L -f + 

z dz 



0, 


(36) 


where Co = x + u. Note that in the new coordinate z, the horizon of the BH is located at 
z — 0 and infinity at z = 1. With the change <p{z) = z a ( 1 — z) l3 F(z), the last equation 
reduces to the hypergeometric differential equation for the function F(z), that is, 


z( 1 - z)F"(z) + (c — (a + b + l)z)F'(z) - abF(z) = 0. (37) 

In this case the coefficients a, b and c are given by the relations 

c = 2a + 1, (38) 

a + b = 2a + 2/3 + 1, (39) 

ab = a (a — 1) + /3(/3 — 1) + 2a + 2/3 + 2 a/3, (40) 

providing the expressions for the coefficients 

a = a + (3, (41) 

b = 1 + cr + /I, (42) 

and for the exponents we obtain 

a = (3 = —ix + u>. (43) 


Without loss of generality, we have chosen the negative signs for the exponents. The solution 
of the radial equation reads 

F(z) = CiFi(a, b, c; z) + C 2 z 1 ~ c F 1 (a — c + 1, b — c + 1, 2 — c; z), (44) 
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where C\ and C 2 are arbitrary constants and F\ (a, b, c; z) is the hypergeometric function. 
The solution for <p(z) is given by 


<p(z) = C lZ - ix+w (l - ^)-“+ w F 1 (a, 6, c; z) + C 2 z ix +“( 1 - z)~ ix+UJ F^a-c+l,b-c+l,2-c]z). 

(45) 

In the neighborhood of the horizon z — 0, the function if(z) behaves as 


tp(z) = Cie-“+“'“ s + C 2 e“+" ln ‘, 


(46) 


for the scalar field 0 one gets 


0 ~ C ie -^ t+x + lnz ^ + Coe~ iu< ' t ~ x+lnz \ 


(47) 


The first term in the last equation corresponds to an ingoing wave at the BH, while the 
second one represents an outgoing wave. In order to compute the QNMs, we must impose 
that there exist only ingoing waves at the horizon of the BH, then C 2 = 0. The radial 
solution at the horizon is given by 

<p(z) = C lZ - ix+UJ (l - z)~ ix+UJ Fi(a, b, c; z). (48) 


In order to implement the boundary conditions at infinity, z — 1, we use the linear trans¬ 
formation z —> 1 — z, and then we apply Hummer’s formula 42] for the hypergeometric 
function, 

<p(z) = C lZ ~ ix ^( 1 - z )-^ ^ c)r(c ~ a ~^ F 1 (a, b; a + b - c + 1; 1 - z) 

1 (c — a)I (c — 0 ) 

+ C lZ - ix +“( 1 - z) ix +“ T ^ a ~t. b t ~ c) F 1 (c - a, c - b- c - a - b + 1; 1 - z), (49) 

T(a)T(6) 

This solution near the infinity, z — 1, takes the form 

<p(z) = Cl (l - 2 )-.^ r(c )r (c-a- 6 ) + Ci(1 _ , (50) 


T(c — a)T(c — b) 
and the scalar held solution near inhnity behaves as 


T(a)T(b) 


0 ~ c L c~ iuj{t+x+ H1 ~ z)) _ a ~ S + ClC -Mt-x+ ind- Z )) r (c) r (q + b- c) 

T(c-a)T(c-6) T(a)T(6) 


To compute the QNMs, we also need to impose the boundary conditions on the solution of 
the radial equation at inhnity, meaning that only purely outgoing waves are allowed there. 
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Therefore, the second term in the last equation must vanish; this is fulfilled, at the poles of 
T(a) or T(6), where the scalar field satisfies the considered boundary condition only when 


a = —n or b = —n, (52) 

where n = 0,1, 2,.... These conditions determine the form of the quasinormal modes, 


CO = 


2x A 


(n+ 1). 


(53) 


B. Spacetime metric ds 2 = 


(x 2 + 2 Mx — e) dt 2 + — - dx 2 

VVV ) (|)i 2 +2Mi-e 


For the second metric given in Eq. (ED, we have a spacetime that is not asymptotically 
flat; as such, and as mentioned before, we use a definition for the QNMs different from 
the one used in an asymptotically flat spacetime. The formal treatment for this kind of 
spacetime is discussed in [3], where they defined QNMs to be modes with only ingoing 
waves near the horizon and vanishing at infinity. Thus, the Klein-Gordon equation (ITUD can 
be written as 


1 2 1 

d t d t 0 + —ud u (j) + 7 y(w 2 - v 2 + )d u d u (j) - m 2 (j) = 0, 


u 2 -uV z ~^ ' 1 2 ~~ u ' r ' l 2 


(54) 


now, we will consider a solution of type (j>(t,u) = ip(u)e lut and definition l = (A/??) 1 / 4 , for 
which the radial equation can be written as follows 


^(w 2 - u\)ip"(u) + ^ uip\u ) + 


or 


u 2 — u\ 


— m 2 I cp(u) = 0, 


(55) 


where the prime denotes derivates with respect the variable u. If we define the change 
of variable z = 1 — u 2 + /u 2 43] and follow the procedure stated for the massless case, the 
equation (1531) transforms into the hypergeometric differential equation (1571) for the function 
F(z), where the coefficients a, b, c are given by the following relations 

1 


ci T b — 2a T 2/3 ~t~ 


2’ 


which gives 


3 3 

ab — a{a — 1) T /?(/? — 1) T 2a(3 T T ~/3, 
c = 2a + 1, 


a = a + /3, 
b = a + (3 + -, 


(56) 

(57) 

(58) 

(59) 

(60) 
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and for the exponents a and ft 


a = —i — uj, (61) 

u + 

P = j(l- Vl + 4 m 2 / 2 ) , (62) 

where, without loss of generality, we have chosen the negative signs. The solution of the 
radial equation reads 

F(z) = CiFi(a, b, c; z) + C 2 z 1 ~ c F 1 (a — c + 1, b — c + 1, 2 — c; z), (63) 


where C\ and C 2 are arbitrary constants and F\ (a, b, c; z) is the hypergeometric function. 
Since (p(z) = z a ( 1 — z)^F(z), the behaviour of the scalar held near the horizon (z = 0) is 
given by 

0 ~ C x e v “+ > + C 2 e ^ u + >. (64) 

Then, the scalar held (f> is purely ingoing at the horizon for C 2 — 0, and therefore the radial 
solution is 

ip(z) = Ciz a (l — z) /3 Fi(a, b, c; z). (65) 

In order to implement boundary conditions at infinity (z = 1), we use the linear trans¬ 
formation z —» 1 — z for the hypergeometric function and we obtain 

(p(z) = C lZ a (l - — J^Fi(a, b] a + b - c + 1; 1 - z) 

1 (c — a)l (c — b) 

+ C lZ a (l - z) c - a - b+ ^ c) ^ a ^. b ~ C) F 1 (c-a,c-b ]C -a-b+l-l-z). (66) 

r (a)T(6) 

Using the condition of the hux 


T ~ tp*(z)d z tp(z) - (p(z)d z (p*(z), 

-2 iC\ — ((1 - z) 1/2 ^~ Vl+4mH2 ) r x 2 + 2 Ti Y 2 ^[YF~ z + (1 - ^ 1 / 2 ( 1 +Vi+4m^p) r2 2 j . 

~ -2 iCl 1 ^- ((1 - z) 2p T, 2 + 2 Ti T 2X /r^I + (1 - z) l ~ 2fi T 2 2 ) , (67) 

where 


T(c)T(c — a — b) 
T(c-a)T(c-6)’ 
T(c)T(a + b — c) 

r(o)r( 6 ) ’ 


( 68 ) 

( 69 ) 
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FIG. 1. In this plot we depict the behaviour of the QNMs expressed in Eq. m, with some 
paremeters values, l = 1 and u+ = y/b. We can see that the BH becomes unstable for large values 
of the mass m. 


then, the flux fl67|) has a leading term (1 — z) 1-2/3 and vanishes at infinity only if we impose 
that 

a = —7i or b = —n, (70) 


where n = 0,1, 2, 3.... These conditions lead to the determination of the quasinormal modes 
as follows 


CO = 



4n + 3 — y/l + 4 m 2 l 2 


(71) 


Our results are represented in Fig. ([TJ) , where it is possible to see that, for a scalar held 
with large mass, the BH becomes unstable, while for lower values of the mass, or m = 0, 
this kind of black hole is stable. 


IV. FINAL REMARKS 


This article was devoted to studying the response of two 1 + 1 BH under scalar pertur¬ 


bations. We focused on the BH solutions found in Ref. 


35] in the context of HL gravity. 
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The BH studied in the present paper also correspond to solutions arising from standard GR 
plus the dilaton held; therefore, the physical properties of these BH can be used in different 
contexts. We noted that in studying the QNM oscillations of the metric (JTTj) with massive 
scalar held perturbations it is necessary to look at the solution in terms of the confluent hy¬ 
pergeometric, or Hummer, functions, and, as a result, found two different cases, in one case 
we have absent QNMs under scalar perturbations and in the second case we have a discrete 


spectrum. These results are different from those obtained in Ref. 


, where the QNMs 


are a continuous spectrum. Also, we computed the frequencies of the massless scalar held 
as sources of perturbations, and again obtained a discrete spectrum. From these results, we 
conclude that this BH is stable under massive and massless scalar perturbations. 

On the other hand, for spacetime in which the cosmological constant does not vanish, we 
found, in addition to the exact quasinormal frequencies, that it is possible to see that when 
the mass of the scalar held is large, the geometry becomes unstable. Finally, we would like 
to note that the frequencies found in this article are purely imaginary, and as such represent 
pure damping behaviour. 


ACKNOWLEDGMENTS 

The authors would like to thank Samuel Lepe and Olivera Miskovic for useful comments. 
MC was supported by PUCV through Proyecto DI Postdoctorado 2015. MG-E acknowledges 
support from a PUCV doctoral scholarship. 


[1] P. Horava, “Quantum Criticality and Yang-Mills Gauge Theory,” Phys. Lett. B 694, 172 

( 2010 ). 

[2] P. Horava, “Membranes at Quantum Criticality,” JHEP 0903, 020 (2009). 

[3] P. Horava, “Spectral Dimension of the Universe in Quantum Gravity at a Lifshitz Point,” 
Phys. Rev. Lett. 102, 161301 (2009). 

[4] P. Horava, “Quantum Gravity at a Lifshitz Point,” Phys. Rev. D 79, 084008 (2009). 

[5] S. Hod, “Bohr’s correspondence principle and the area spectrum of quantum black holes,” 
Phys. Rev. Lett. 81, 4293 (1998). 


15 




[ 6 ] K. D. Kokkotas and B. G. Schmidt, “Quasi-normal modes of stars and black holes,” Living 
Rev. Rel. 2, 2 (1999). 

[7] G. T. Horowitz and V. E. Hubeny, “Quasinormal modes of AdS black holes and the approach 
to thermal equilibrium,” Phys. Rev. D 62, 024027 (2000). 

[ 8 ] V. Cardoso and J. P. S. Lernos, “Quasi-normal modes of Schwarzschild anti-de Sitter black 
holes: Electromagnetic and gravitational perturbations,” Phys. Rev. D 64, 084017 (2001). 

[9] V. Cardoso, R. Konoplya and J. P. S. Lemos, “Quasinormal frequencies of Schwarzschild black 
holes in anti-de Sitter spacetimes: A complete study on the asymptotic behavior,” Phys. Rev. 
D 68 , 044024 (2003). 

[10] J. Natario and R. Schiappa, “On the classification of asymptotic quasinormal frequencies for 
d-dimensional black holes and quantum gravity,” Adv. Theor. Math. Phys. 8 , 1001 (2004). 

[11] V. Cardoso, J. Natario and R. Schiappa, “Asymptotic quasinormal frequencies for black holes 
in non-asymptotically flat spacetimes,” J. Math. Phys. 45, 4698 (2004). 

[12] J. S. F. Chan and R. B. Mann, “Scalar wave falloff in topological black hole backgrounds,” 
Phys. Rev. D 59, 064025 (1999). 

[13] B. Wang, E. Abdalla and R. B. Mann, “Scalar wave propagation in topological black hole 
backgrounds,” Phys. Rev. D 65, 084006 (2002). 

[14] R. A. Konoplya, “Gravitational quasinormal radiation of higher-dimensional black holes,” 
Phys. Rev. D 68 , 124017 (2003). 

[15] P. Kanti and R. A. Konoplya, “Quasi-normal modes of brane-localised standard model fields,” 
Phys. Rev. D 73, 044002 (2006). 

[16] J. S. F. Chan and R. B. Mann, “Scalar wave falloff in asymptotically anti-de Sitter back¬ 
grounds,” Phys. Rev. D 55, 7546 (1997). 

[17] V. Cardoso and J. P. S. Lemos, “Scalar, electromagnetic and Weyl perturbations of BTZ black 
holes: Quasi normal modes,” Phys. Rev. D 63, 124015 (2001). 

[18] J. Crisostomo, S. Lepe and J. Saavedra, “Quasinormal modes of extremal BTZ black hole,” 
Class. Quant. Grav. 21, 2801 (2004). 

[19] E. Berti, V. Cardoso and J. P. S. Lemos, “Quasinormal modes and classical wave propagation 
in analogue black holes,” Phys. Rev. D 70, 124006 (2004). 

[20] S. Lepe and J. Saavedra,“Quasinormal modes, superradiance and area spectrum for 2+1 
acoustic black holes,” Phys. Lett. B 617, 174 (2005). 


16 



[21] J. Saavedra,“Quasinormal modes of Unruh’s acoustic black hole,” Mod. Phys. Lett. A 21, 
1601 (2006). 

[22] V. Ferrari, M. Pauri and F. Piazza, “Quasi-normal modes of charged, dilaton black holes,” 
Phys. Rev. D 63, 064009 (2001). 

[23] S. Fernando and K. Arnold, “Scalar perturbations of charged dilaton black holes,” Gen. Rel. 
Grav. 36, 1805 (2004). 

[24] R. A. Konoplya, “Decay of charged scalar field around a black hole: Quasinormal modes of 
R-N, R-N-AdS and dilaton black hole,” Phys. Rev. D 66, 084007 (2002). 

[25] S. P. Robinson and F. Wilczek, “A relationship between Hawking radiation and gravitational 
anomalies,” Phys. Rev. Lett. 95, 011303 (2005). 

[26] Y. S. Myung and H. W. Lee, “Schwarzschild Black Hole In The Dilatonic Domain Wall,” 
Phys. Rev. D 63, 064034 (2001). 

[27] T. Torii and K. I. Maeda, “Stability of a dilatonic black hole with a Gauss-Bonnet term,” 
Phys. Rev. D 58, 084004 (1998). 

[28] E. Teo, “Statistical entropy of charged two-dimensional black holes,” Phys. Lett. B 430, 57 
(1998). 

[29] M. D. McGuigan, C. R. Nappi and S. A. Yost, “Charged black holes in two-dimensional string 
theory,” Nucl. Phys. B 375, 421 (1992). 

[30] R. Becar, S. Lepe and J. Saavedra, “Quasinormal modes and stability criterion of dilatonic 
black hole in 1+1 and 4+1 dimensions,” Phys. Rev. D 75, 084021 (2007). 

[31] R. Becar, S. Lepe and J. Saavedra, “Decay of Dirac fields in the backgrounds of dilatonic 
black holes,” Int. J. Mod. Phys. A 25, 1713 (2010). 

[32] A. Lopez-Ortega and I. Vega-Acevedo, “Quasinormal frequencies of asymptotically flat two- 
dimensional black holes,” Gen. Rel. Grav. 43, 2631 (2011). 

[33] R. Becar, P. A. Gonzalez and Y. Vasquez, “Dirac quasinormal modes of two-dimensional 
charged Dilatonic Black Holes,” Eur. Phys. J. C 74, 2940 (2014). 

[34] D. Bias, O. Pujolas and S. Sibiryakov, “A healthy extension of Horava gravity,” Phys. Rev. 
Lett. 104, 181302 (2010). 

[35] D. Bazeia, F. A. Brito and F. G. Costa, “Two Dimensional Horava-Lifshitz Black Hole Solu¬ 
tions,” Phys. Rev. D 91, no. 4, 044026 (2015). 


17 



[36] R. B. Mann, S. M. Morsink, A. E. Sikkema and T. G. Steele, “Semiclassical gravity in (1+1)- 
dimensions,” Phys. Rev. D 43, 3948 (1991). 

[37] D. Christensen and R. B. Mann, “The Causal structure of two-dimensional space-times,” 
Class. Quant. Grav. 9, 1769 (1992). 

[38] S. Estrada-Jimenez, J. R. Gomez-Diaz and A. Lopez-Ortega, “Quasinormal modes of a two- 
dimensional black hole,” Gen. Rel. Grav. 45, 2239 (2013). 

[39] James B. Seaborn, “Hypergeometric Functions and Their Applications,” Springer-Verlag, New 
York (1991). 

[40] A. Farrell and B.P. van Zyl, “Universality of the energy spectrum for two interacting harmon¬ 
ically trapped ultra-cold atoms in one and two dimensions,” J. Phys. A: Math. Theor. 43, 
015302 (2010). 

[41] A. D. MacDonald, “Properties of the confluent hypergeometric function,” Technical Report 
No. 84, Research Laboratory of Electronics, Massachusetts Institute of Technology, 1948. 

[42] M. Abramowitz and A. Stegun, Handbook of mathematical functions, (Dover Publications, 
New York, 1970). 

[43] R. Aros, C. Martinez, R. Troncoso and J. Zanelli, “Quasinormal modes for massless topological 
black holes,” Phys. Rev. D 67, 044014 (2003). 


18 



